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EFFECTS OF LEADING-EDGE BLUNTNESS ON 


PRESSURE AND HEAT-TRANSFER MEASUREMENTS OVER A 

FLAT PLATE AT A MACH NUMBER OF 20 

By William D. Harvey 
Langley Research Center 


SUMMARY 


An experimental investigation has been conducted to study the effects of 
viscous interaction and leading-edge bluntness on flat-plate pressures and heat- 
transfer coefficients at a Mach number of approximately 20. The investigation 
was conducted at a nominal stagnation temperature and pressure of about 4000° R 
and 10 000 pounds per square inch for a free- stream Reynolds number per foot of 
about 0.b6 x 10^. 

The pressure data for the sharp leading-edge flat plate indicated a region 
of transition between the free molecule flow and continuum flow regions near 
the leading edge. Experimental results from the blunt leading-edge configura- 
tions were In agreement with those from the linear addition of inviscid and 
viscous effects. The blast-wave theory gave good approximations of the shock 
shapes of the blunt leading-edge configuration. The boundary-layer growth meas- 
ured from schlieren photographs of the sharp plate gave fair agreement when 
compared with the hypersonic boundary- layer theory. 


INTRODUCTION 


There is continuing interest in low-density hypersonic flow over simple 
lifting shapes, such as flat plates, to study the effects of leading-edge blunt- 
ness and viscous-layer thickness on the distributions of pressure and heat- 
transfer rates over the plates. A number of theoretical and experimental inves- 
tigations have been made to study these effects. Most of these investigations 
consider extremely blunt leading edges for bluntness effects or extremely sharp 
leading edges for boundary-layer — shock-wave interaction effects. The blast- 
wave theory has been used to predict the dominant bluntness effects with a fair 
amount of success (refs. 1 to 6) , although a number of shortcomings of the 
blast-wave analogy are shown in references 7 and 8. The well-established weak 
and strong interaction theories have been used extensively to predict the 
continuum-flow region of the do m inant boundary- layer displacement effects on 
the sharp leading-edge flat plate (refs. 9 to 21). At the present time, no 
theory adequately predicts the transition region between continuum flow and 



free molecule flow near the leading edge of the plate ; however, some data have 
been presented which include slip flow effects. (See for instance, refs. 22 
to 2k). 

Several experimental investigations have been made to study hypersonic 
flow over flat plates with blunt leading edges (refs. 25 to 50 ) . Some of these 
studies were conducted with adiabatic wall conditions for which the boundary- 
layer displacement effects are considerably larger than those anticipated for 
cold- wall tests considered in hypersonic -flow phenomena (refs. 25 and 26 ). 

Since the investigations often did not include both pressure and heat-transfer 
measurements, the analysis of the experimental data was limited. 

The purposes of the present investigation were to measure the surface- 
pressure distribution and the heat -transfer-rate distribution simultaneously 
and to obtain shock-wave shapes from schlieren pictures for low-density hyper- 
sonic flow over sharp and blunt leading -edge configurations. This investigation 
was made in the Langley hotshot tunnel with nitrogen at a nominal free -stream 

Mach number of 20, a nominal Reynolds number per foot of 0.k6 x 10^, and a 
stagnation temperature of 1000° R; the corresponding ratios of wall temperature 
to stagnation temperature were approximately 0.13- The plates were alined with 
the stream at angles of attack and yaw of 0°. The variation in leading-edge 
bluntness from 0.0003 inch to 0.080 inch yielded Reynolds numbers based on 
leading-edge thickness from approximately 11. 5 to 3060 . 

Shock shapes were determined from schlieren photographs, and detailed sur- 
face heat-transfer and pressure measurements were made. Correlation of the 
experimental data was made in terms of local similarity. Also, a comparison of 
the extended viscous similitude for bluntness was made along with a discussion 
of data and theory agreement. 


SYMBOLS 


A, B,D coefficients in equations for the zero-order strong-interaction 
theory (from ref. 16) 

Cp nose drag coefficient 

C w viscosity-temperature ratio, tV^coA^w 

C' viscosity-temperature ratio based on T* method, lAo^oo/^oo^-oo 

c specific heat of calorimeter sensing disk material, Btu/lb-°F 

Cp w free-stream specific heat at constant pressure 

E electromotive force, volts 

0 = f(T w ,T r ,7,Mb)(from ref. 21 ) 
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local heat-transfer coefficient without pressure gradient, 

(dp/dx = 0) 

local heat-transfer coefficient with pressure gradient included 
(dp/dx 4 0) 

coefficients defined by reference 16 

thermocouple gage -calibration constant, millivolts/°F 
free-stream Mach number 

mass, pounds 

Prandtl number at free-stream static temperature 

Stanton number including pressure gradient effect, based on free 

stream properties, - 

C p,coR»X» 

exponent in equation for variation of pressure with x, p « x n , 
the slope of the power-law curve 

pressure, lb/sq in. 

wall pressure, lb/sq in. 

free-stream pressure, lb/sq in. 
heat content, Btu/sq ft 

local heat-transfer rate, dq/dT, Btu/ft^-sec 
free-stream Reynolds number per foot 

free-stream Reynolds number with x as characteristic length, 
Pco V oo X 
Poo 

Reynolds number based on leading-edge thickness 

radius of calorimeter disk 
area, sq in. 
temperature, °R 

calculated arc chamber temperature, °R 


also 


3 



T 


T, 


w 


Tl 


t 

V 

V„o 

X 

y 


z 5 

a 

7 


P 

P 

arc 


T 

e 

♦ 


recovery temperature, °R 


wall temperature, °R 

stagnation temperature behind normal shock, °R 
free -stream static temperature, °R 
temperature evaluated by T* method (ref. 31) 
thickness, inch 
volume, cu in. 

free-stream velocity, ft/sec 
distance along plate surface 

spanwise distance on plate from model center line 

vertical distance from origin (center line through model leading-edge 
thickness, t^g) to shock wave (See sketch A) 

vertical distance from plate surface to outer edge of visible boundary 
layer (See sketch A) 

angle of attack, degrees 

ratio of specific heats 


hypersonic viscous parameter. 



sensing -material density, lb/ft3 
arc chamber density, atm 


free-stream density, lb/ft 5 


free-stream viscosity; viscosity at the wall 
time, sec 

correction factor in equation (5) 
angle of yaw, deg 
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Subscripts : 

c calorimeter disk 

Ze leading edge 


APPARATUS AND TESTS 


Tunnel and Test Conditions 

The present investigation was conducted in the Langley hotshot tunnel. 

This tunnel is similar to other hotshot facilities in that a high-energy arc is 
discharged within an arc chamber to heat and pressurize the test gas. Upon rup- 
ture of a diaphragm upstream of the nozzle throat, the nitrogen expands through 
a conical nozzle and test section into the vacuum reservoir. A detailed 
description of the tunnel can be found in reference 32. The arc chamber used 
in the present tests is the coaxial electrode arrangement described in appen- 
dix A of reference 32 without the magnetic coil. The present investigation was 
made at a nominal Mach number of 20, and arc chamber temperatures and pressures 
of approximately k000° R and 10 000 pounds per square inch with nitrogen as the 
test medium. 

During each test, the arc chamber pressure and test-section pitot pressure 
were measured and recorded. Two transducers for each measurement were employed 
and relatively good agreement percent) existed between the two transducers 
for each case. Calibrations in the anticipated pressure range for the particu- 
lar test were made on each transducer beforehand. 


Schlieren Apparatus 

The schlieren system employed was a conventional horizontal off-axis 
single-pass type with 8 -foot -focal -length 12-inch-diameter parabolic mirrors. 

The light source was a mercury arc lamp designed for either continuous or short- 
duration arc service. The short -duration arc light was programed at a prese- 
lected time for each test of the present investigation. Schlieren photographs 
were obtained on 4--inch by ^-Inch. film plates. Model roll and yaw were care- 
fully adjusted to match the upper rays in the schlieren beam. The angle of 
attack was set with an inclinometer and was equal to zero for all 
configurations . 


Model 

The flat-plate model employed in the investigation is shown in figure l(a). 
For the tunnel tests, the screw holes were filled. Figure l(b) is a photograph 
of the underside of the instrumented surface and shows the pressure orifices 
and reference manifold. The pressure orifice tubes are shown without pressure 
transducers. Thermocouples are in place. The large cylinder with small tubes 
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protruding outwardly is the reference pressure manifold common for all 
transducers . 

Figure 2 is a schematic drawing of the model and shows the pertinent 
details of the model construction and instrumented surface. Table I gives the 
location of the instrumentation. The model housing and sting are made of alumi- 
num except for the instrumented plate surface which is made of stainless steel 
and is discussed subsequently in more detail. There are four leading -edge 
pieces semieylindrical in shape and interchangeable. Leading-edge thicknesses 
are 0.000 3 inch, 0.015 inch, 0.030 inch, and 0.080 inch. Each leading-edge 
piece was measured with a microscope and the average thickness across the span 
of the model was determined with a ±0. 0001-inch variation. 


Heat -Transfer Gages 

Heat-transfer rates of 1 to 10 Btu/ft^-sec were expected during the present 
investigation with tunnel test times of about 0.100 second. In order to develop 
a gage that would record a nearly instantaneous surface -temperature change with 
time and within a very small error range, it was necessary to consider the 
following equation that expresses the quantity of heat contained by any body: 


<1 


mcT 

S 


( 1 ) 


The rate of heat transfer to the body is given by 


4 


dq 

dT 


( 2 ) 


Therefore, by differentiation and substitution of equation (l) into equation (2) 


• me dT 

q S dr 


(3) 


where m = pV. For a disk V = and S = nr^. Substituting these equiv- 

alents into equation (3) yields 
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me dT _ pVc dT _ dT 

S dr ~ S dr ” pC c dT 


00 


The change in the electromotive force in a thermocouple may be expressed 
as 


dE = k dT 

where k is a constant equal to the change in millivolts per degree Fahrenheit 
for a particular thermocouple. By substitution of dE = k dT into equa- 
tion (4), the following expression is obtained 
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Since k, p, and c are essentially constant for a particular choice of mate- 
rial, only the thickness t c need be varied to adjust the sensitivity. Thus 
by reducing the thickness for a given thermoelement, the sensitivity is 
increased. 

In constructing a calorimeter thermocouple gage for making the desired 
heating-rate measurements, the following factors were considered: 

(1) The calorimeter was constructed so that the anticipated temperature 
rise would not exceed that of the working limits of the sensing material, 
bonding material, or support structure. 

(2) The ratio of sensing-disk radius to thickness was made large to pre- 
vent heat losses radially outward from the material sensing area. 

( 3 ) The calorimeter needed sufficient sensitivity to produce an electrical 
signal which could be recorded with accuracy. 

The heat-transfer gage unit was constructed at the Langley Research Center. 
A support plate was made of steel and was initially machined and ground on one 
side. Then 20 holes for pressure-orif ice tubes (0.060-inch -diameter) and 
21 0 . 25 - inch- diameter holes for thermocouples were drilled at specified loca- 
tions. A 0.001- inch-thick bonding material similar to double-backed tape was 
used to attach the 0.002-inch-thick 302 stainless-steel sheet to the support 
plate. This bonding resulted in a smooth uniform thin-skin surface, which 
served as a calorimeter-type heat-transfer gage at each of the 0 . 25 -inch hole 
locations. Chromel-alumel thermocouple wires (0. 001-inch diameter) were then 
resistance-welded on to the back side of the thin skin as near the center of 
the disks as possible. The wires were limited to l/4 inch in length to minimize 
resistance. Copper lead wires of much larger diameter were joined to the 
0.001-inch-diameter thermocouple wire which served as the cold junction for 
each gage. 

Calibration of the gages was made by applying a known electromotive force 
to each thermocouple circuit including the readout equipment. The variation of 
the deflection in inches with millivolts was recorded for each gage prior to a 
test. Since the sensitivity was greatly increased by reducing the calorimeter 
thickness, it was found that with the proper galvanometer selection, no ampli- 
fication was necessary to record the thermocouple output. The repeatability of 
the electromotive force calibration was about 2 percent for the linear output. 


Model Pressure Gages 

Pressure measurements on the flat plate were made with single -diaphragm 
variable -reluctance -type transducers. They were rated at 0.25 pound per square 
inch full scale with an excitation voltage of 5 volts at 20 kilocycles. The 
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response time of the transducers during the present investigation was well 
within the order of 1 millisecond at full scale as specified by the 
manufacturer . 

The transducers were mounted in the support plate so that the orifice tube 
length was minimized to 3/8 inch to avoid unnecessary lag in the measurements . 
(See fig. 2.) Carrier amplification of the transducer output signal allowed 
pressure measurements as low as 0.005 psi to be made with reasonable accuracy. 
Calibration of each model pressure transducer was made over the pressure range 
anticipated for the test. The linear calibrations could be repeated within 
±5 percent. 

Table I shows two pressure orifice locations on each respective leading- 
edge piece (not shown in fig. 2). One pressure orifice was located at the same 
x- and y-coordinates on each leading -edge piece. The second orifice had the 
same y-coordinate but different x-coordinates for each leading edge. 


DATA REDUCTION 


Heat -Transfer Measurements 

In order to obtain, reduce, and analyze the heat -transfer data, it was 
necessary to make the following assumptions, some of which are based on pre- 
vious experience: 

(1) The sensing skin was thin enough that the temperature of the under- 
surface was representative of the average temperature through the skin thick- 
ness. Calculations for the 0.002-inch-thick stainless -steel sensing surface at 
a stagnation point on a model indicated that the back-surface temperature was 
equal to that of the average front temperature, a thermal diffusivity of 

5 x 10*5 ft 2 /sec being considered. 

(2) Heat losses from the periphery of the calorimeter sensing surface to 
the support structure were assumed to be negligible. 

(3) Heat applied to the sensor entered the front surface only. 

(k) Uniform heat distribution existed over the entire calorimeter surface. 

(5) Constant thermal diffusivity existed throughout the sensing surface. 

(6) The thermocouple wires sensed the back-surface temperature of the 
material . 

(7) Heat loss along the length of the thermocouple wires was negligible. 
Calculations were made by assuming that the two wires attached to the underside 
of the sensing material were of equal diameter. These calculations indicated 
that for a constant q there was a heating-rate error for the two wires of 
about 3 percent at 20 milliseconds after heat was initially applied to the 
sensing surface. 
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(8) No heat sink existed at the junction of the thermocouple wires and 
sensing underside surface. 

The factors which led to the validity of the preceding assumptions were 
as follows: 


(1) The thin stainless-steel skin was precision rolled to a thickness of 
0.002 inch ± 0.00005 inch. 

(2) Handbook values of p and c were used, and the temperature rise of 
the disk material was kept to a few degrees above room temperature. Therefore, 
the density and specific heat were essentially constant. 

( 5 ) A linear thermoelement output curve was used. 

(4) Since the diameter (0.25 inch) of the sensing disk was much larger 
than the thickness ( 0.002 inch), it was felt that the heat loss radially out- 
ward or to a sink at the disk periphery was negligible. 

( 5 ) The temperature gradient through the sensing surface was negligible 
and no loss due to wire length was experienced. It is believed that the heat- 
transfer slopes read from the records were accurate within 1 percent. 


Therefore, since the thickness t c (0.002 inch), specific heat c 

(0.12 Btu/lb-°F), and density of the sensing material ( 0.29 lb/in5) were known, 
it may be seen from equation (^) that only the time derivative of the average 
surface temperature dT/dr was necessary to obtain heat-transfer data. This 
derivative was actually determined by measuring the slope of each thermocouple 
output trace at discrete time intervals in terms of inches of deflection. The 
paper speed of each oscillograph record (ft/sec) together with the change in 
voltage per degree Fahrenheit for the chromel-alumel thermocouple wires 

( — to — E and the calibration slope for each gage circuit 


\ millivolts/ 

used to compute general heat-transfer-rate equation, as follows 


ft 


were 


q = pct 0 — = pet (4^) (paper speed) (calibration slope) (trace slope) 
1 dT G 


T>4- -1 ■j 

where 4 is measured in — — . 

ft 2 -sec 


Pressure Measurements 

Prior to testing, a calibration of the pressure range anticipated for the 
test was made for the individual gages. Measurements were read from oscillo- 
graph records at 10-millisecond intervals with zero time selected as the first 
indication of a pressure rise. 
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Correction to Data 


Until recently, it had been the practice of investigators in hotshot tun- 
nels to define the flow by simultaneously measuring the test-section pitot 
pressure and the reservoir pressure, and from the known density of the initial 
charge of gas in the reservoir, to assume that the same density, uniformly dis- 
tributed, existed in the reservoir after arc discharge. With these three param- 
eters it was possible to calculate all the other characteristics of the flow 
with a suitable data-reduction program such as that used in reference 33. 

Recent work (refs. 3^ and 35 a 11 *! unpublished data from other facilities) has 
indicated that the assumption of uniform density in the reservoir may lead to 
serious error in flow definition. To define the flow more accurately, a series 
of tests was made for which simultaneous measurements of free-stream velocity, 
stagnation-point heat-transfer rate, and the previously mentioned pressures and 
density were obtained. These results yield correction factors that may be 
applied to test results for which the density was assumed to be uniformly 
distributed. 

Several simultaneous measurements of the free-stream velocity and 
stagnation -point heat -transfer rate were made with the flat plate installed in 
the tunnel. All the data of the present investigation have been corrected by 
the factors which are presented in table II. 


THEORY 


Hypersonic Self-Induced Pressure Gradient on a Sharp 
Leading -Edge Flat Plate 


Several investigations have shown that an asymptotic solution for large 
boundary-layer-induced pressures can be related to a physical plane pressure 
gradient p « x n if hypersonic and isentropic flow is assumed at the boundary- 
layer edge. (See refs. 36 to 3 9 -) 

Bertram and Feller (ref. 16) have presented an analysis which follows that 
of Li and Nagamatsu (ref. 36 ) for laminar boundary- layer flows in which an axial 
pressure gradient exists along the plate surface. Equations are presented in 
reference 16 which show that the induced pressure ratio, boundary-layer thick- 
ness, and local heat-transfer characteristics can be written as: 


Pw _ A *gjjV 
P ~ ^ 


(5) 


fs 

X 




( 6 ) 


10 




where A , B, and D are coefficients given by equations (21), (22), and (2^) 
of reference 1 6 . These coefficients are plotted in figure 9 of reference 1 6 as 
a function of wall-to-stagnati on -temperature ratio for various values of 7 . 
Equations (5) to (7) are applicable for a hypersonic self -induced pressure gra- 
dient on a flat plate for a strong interaction solution. 

Bertram and Feller have also shown in reference 1 6 that from the analysis 
of Li and Nagamatsu (ref. 36 ) for hypersonic flow, where the Prandtl number is 
equal to unity and n (exponent of p « x n ) may have positive or negative 
values, the heat-transfer coefficients may be expressed as the ratio 

(8) 

where h represents the case in which dp/dx = 0. The coefficient is 

defined by equation (9) of reference 16 and accounts for the effect of the 
pressure gradient. The amount that Kj deviates from unity is a measure of 
the importance of this gradient apart from the effect of changes in local con- 
ditions. It was found earlier that was relatively insensitive to changes 

in wall temperature and 7 for the insulated-wall case. The square root of the 
pressure ratio accounts for differences between local and undisturbed free- 
stream conditions. The following sketches are typical of the model used for 
the numerical computations of the preceding theory for boundary- layer displace- 
ment and blast -wave analysis. 
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Pressure Distribution for Blunt Leading-Edge 
Flat Plate at a = 0° 

Inviscid effect .- The inviscid effect of blunt leading edges usually dom- 
inates the induced pressure and heat transfer for high Mach number low-density 
flow for several nose diameters downstream of the leading edge. This effect 
results from an entropy jump at the leading edge near the plate surface. In 
particular, Bertram and Blackstock (ref. 21 ) have shown that the addition of 
the pressure increment obtained from tangent -wedge theory to the increment due 
to blunting the leading edge on a flat plate, gave results which were in good 
agreement with the characteristic theory. (See ref. 6.) Therefore, for the 
case where P w /p 0O = 1.2, the blunt -leading-edge inviscid effects can be 

expressed in the form 


•w _ 


= 0 . 187 e 


fr(r - i) M “ 


X 


/tie] 


2/3 


+ 0.74 


(9) 


where 

e « 1 - [0. 0048(7 - 1)0 

and represents a small correction term to the function of 7 in equation (9)* 

Small viscous effects . - If the blunt -leading-edge induced pressures are 
assumed to be much larger than the viscous induced-pressure increment and some- 
what larger than the free-stream pressure, the tangent-wedge portion of the 
solution can be obtained, as in reference 21, through proper substitution of 
the equation representing the slope of the boundary-layer edge in hypersonic 
similarity form into the first-order equation representing hypersonic small- 
perturbation theory. The completed expression for this solution is 

APviscous _ 5 ^^00 , > 

Poo = 6 U0; 

where for the Prandtl number of O.725 of this investigation 


+ 0.352 j 

If equation (9) is added to equation (l0), an account of the inviscid 
plus weak viscous interaction effects can be evaluated empirically with some 
degree of accuracy. 

Strong viscous effects .- If the viscous -induced increment in pressure is 
too large to be ignored and contributes sufficiently to the induced effect on 
the boundary-layer growth, the effect must be accounted for. Bertram and 


G = 1.6482-=— if 2k 
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Blackstock (ref- 21) have developed an equation which when used with the invis- 
cid case should predict with some degree of accuracy the effects of inviscid 
and strong viscous interaction. By the substitution of the strong interaction 
equation in terms of boundary -layer slope (eq. (15) of ref. 21 ) into the expres- 
sion for the hypersonic similarity deflection angle (based on the edge of the 
boundary layer ), the equation for strong viscous effects becomes 


^viscous 


1 

2 



iii ± 1 ) 
2 





(H) 


Therefore, to account for both the inviscid and strong viscous effects, equa- 
tion (ll) may be added to equation ($))• In equation (ll), n is the exponent 
in the power equation for the surface pressure variation with x, P w « x n . 

In the equation of the hypersonic-similarity theory is a coefficient 

varying with laminar boundary -layer growth, and is expressed as a function of 
pressure gradient and wall-to-stagnation temperature ratio at various values 
of 7 . (See fig. 1 of ref. 21.) 


Shock-Wave Shapes 

Many approaches to the blast-wave concept have been made (refs. 1 to 6 ). 

For inviscid hypersonic flow over a simple slender two-dimensional body, 
Lukasiewicz (ref. l) has shown that approximate solutions of blast analogy can 
be derived for plane or axisymmetric flow. Reference 1 has shown that for 
7 = 1.4 and plane flow at zero angle of attack, the following expression exists 
for shock-wave shapes to the first approximation 
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Correlation of Experimental Data 


A theoretical correlation value for the undisturbed flat-plate Stanton 
number may be obtained by various means where dp/dx = 0. (See refs. 1 6 , 25, 
and 38 . ) In this presentation, the T* method and modified Reynolds number 
analogy were applied to the Blasius skin-friction value to obtain the correla- 
tion parameter for the Stanton number (ref. ^0). 

- 0- 352Np r; ; 2 / 5 R^ x - l/ 2 (C' ) l/a (Ilf) 




where C 1 = — °° and where „ 

u T > Pr,oo 

^00 OO 

evaluated at T' . (See ref. 31*) 


-2/3 


is the Reynolds number analogy factor 


To evaluate the viscosity-temperature ratio, Monaghan expresses the T 
equation in reference 31 as 


mi iji 

-= = Jd. + 0.468U- 

T. Pr >“ V 


1/3 (fe ' fe) ' °- 275H pr,-( z i J; ) M “ (15> 


Bertram and Feller (ref. 16) have further shown that for very high Mach 
numbers (Np r „ = Constant), equation ( 15 ) becomes 

~ = ^(l - 0.468N Pr ^ 1 / 3 ) + 0 A 68 N Pr ^ 5/6 - 0 . 273 N Pr ^ ( 16 ) 

*fc "t 

From equation (l 6 ), the T 1 temperature becomes constant for a given 
ratio of wall temperature to stagnation temperature. In this presentation, a 
value of O .72 was used for the Prandtl number and equation (l 6 ) becomes 

= 0.160 + 0 . 580 ^ ( 17 ) 
T t T t 


If the assumption is made that equation (l^-) is valid when local conditions 
are utilized along with constant wall pressure, the blunt flat plates can be 
simply evaluated by assuming that the flow over the plate has passed through a 
normal shock and the flow at the boundary-layer edge is at constant pressure 
(the sharp leading -edge plate is considered as having a degree of bluntness). 

To account for the case where dp/dx / 0, the heat-transfer-coefficient ratio 
which has been given in equation (8) as 
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is used for hypersonic flow with a Frandtl number of 1 (Li and Nagamatsu; 
ref, 36 ), for positive or negative values of n, and h representing the value 
for which the pressure gradient is zero (dp/cbc = 0), To evaluate the sharp or 
blunt leading-edge case, equation ( 8 ) was multiplied by equation (l^) to obtain 


%t^oo\po^x 




2 /5 5 p o 

"Pr, 00 


(18) 


The measured pressure distributions in this expression can also he used to give 
an independent check of the measured and calculated heat-transfer distributions. 
Equation (l 8 ) may he used to determine an error existing in the heat-transfer 
parameter. In each case, the experimental pressure data were fitted according 
to the power law, p x n . 


RESULTS AND DISCUSSION 


General 

Figure 3 shows the typical variation and reproducibility of the measured 
local heat-transfer and pressure distribution for two successive runs on the 
sharp flat plate at angles of attack and yaw of 0°. In general, the accuracy 
and repeatability of the data were satisfactory. 

At least two tunnel tests were made for each of the leading-edge blunt- 
nesses. The magnitude and distribution of the heat transfer indicate that the 
boundary layer was laminar. Because of somewhat higher heat-transfer rates at 
the most downstream stations it is suspected that in some instances there may 
have been some type of trailing -edge influence on the last few gage locations. 
The discussion regarding spanwise variations in pressure and heat transfer is 
very brief since the small model volume allowed only a limited number of gages 
to be installed. Therefore, the two streamwise rows (see fig. 2) are. given 
emphasis . 

A nominal set of flow conditions for theory was used for the test condi- 
tions noted on the figures since the conditions tested were nearly the same. 
Thus, there exists in each case only one theoretical curve rather than two. 

Spark schlieren photographs were taken during each test. In each case, a 
no-flow picture was taken to determine flaws in the test-section schlieren win- 
dows. The windows produced negligible image distortion. In some cases mis- 
alinement of the schlieren system was suspected, so that the true shock standoff 
distance could not be determined. Typical schlieren photographs are shown in 
figure k. Flaws in the window are evident in the photographs. The outer edge 
of the boundary layer can be distinguished in figure ^(a) for the sharp leading 
edge (t = 0.0003 inch). As the leading-edge bluntness is increased, the density 
decreases near the plate surface. With sufficient blunting, the boundary layer 
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■becomes undetectable as shown in figure 4(d). The apparent reduction in schlie- 
ren sensitivity contributed to the decrease in the ability to detect the bound- 
ary layer in figure 4. 

Sharp flat plate . - Figure 5 shows the comparison between theory and exper- 
iment in terms of the pressure ratio (at top of figure) and the heat-transfer- 

correlation parameter %t,«4 R oojX Tor Two for the sharp flat plate at 

zero angle of attack. The variations of the pressure ratio and the heat- 
transfer -correlation parameter with Reynolds number based on the distance from 
the leading edge along the plate surface are plotted in figure 5* The open 
symbols represent a spanwise location of y = ±0 .75 inch from the center line 
and the solid symbols are for ±0.25 inch from the center line. (See fig. 2 
and table I.) Pressure data for gages 1 6 and 17 (for x = 5-50 inches and 
y = 1.25 inches and y = 1.75 inches) were not plotted and neither were the 
data for heat-transfer gages l6, 17, and l8 (for x = 5-50 inches and 
y = 1.25, y = 1.75* &nd y = 2.25 inches). The data, in general, showed a 
spanwise decrease at these stations. There seems to be no spanwise effect 
within the range represented by the open symbols. 

The induced pressure theory (eq. (5)) was applied and was found to agree 
with the data over most of the sharp plate, although it was somewhat higher 
near the leading edge. The discrepancy between theoretical and experimental 
results near the leading edge is believed to be due to a transition between con- 
tinuum and free -molecular flow since the Knudsen number for the present experi- 
mental investigation on the sharp leading-edge plate was calculated to be about 
2.5. Some other tests were made on the sharp plate in addition to those pre- 
sented herein to give some insight into the extent the transition flow existed 
on the sharp leading-edge plate. With an increase in Knudsen number the 
reduction in pressure and heat transfer near the leading edge of the plate was 
experienced over a longer length of the plate and an increase in the slip or 
transition region was indicated. These preliminary results are similar to those 
of references 22 to 24 and 26. A slope of n = -0.500 was used for the strong- 
interaction theory. A value for in equation (8) was obtained for 7 = 7/5 

from figure 6 of reference l 6 . 

The strong-interaction theory for the pressure distribution (which is in 
good agreement with the pressure data) was used in equation (l 8 ) to calculate 

Ngt,oo^x>,x and - this strong-interaction theory is compared with the experi- 
mental data in the lower plot of figure 5* Both the experimental data and the 
strong-interaction theory have a somewhat higher trend than the hypersonic 
similarity theory (eq. ( 7 ))^ which is also shown in figure 5. This result is 
not surprising since equation (l8) is a more general form of equation (7). 

(See ref. l 6 .) The results from sharp leading-edge theory for the Stanton num- 
ber correlating parameter, given by equation (l4) ^ zero-pressure-gradient 

theory, = p^, ^ = 0 ), were also compared with the experimental data and 

illustrate the effect of pressure gradient on hypersonic heat transfer. 
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Blunt flat plate . - Equation (l8) was used to correlate the experimental 
data plotted in figures 6, J, and 8. Again, the pressure ratio, plotted at the 
top of each figure, was fitted by the power law and the values of n were used 
to obtain Kj. 

Estimations of the experimental pressure distribution for each leading- 
edge bluntness were obtained by individually adding the results from equa- 
tions (10) and (ll) to those of equation ( 9 ). For each bluntness, the weak 
viscous effects (eqs. ( 9 ) and ( 10 )) fell below the experimental data near the 
leading edge, but approached good agreement at the most downstream stations for 
increased leading -edge bluntness. In all cases, the strong viscous effects 
(eq. ( 9 ) and eq. (ll)) agreed well with the experimental pressure distribution, 
except near the leading edge where the theoretical result is lower than experi- 
mental results predicted. Experimental pressure data in reference 26 were 
obtained for test conditions similar to those of the present investigation for 
Mach numbers of 18.3 and 18.7 and 110 < Rj. < 63 O; however, the ratio of wall 
temperature to stagnation temperature was somewhat smaller and the Knudsen num- 
ber larger in reference 26 . The data of reference 26 showed good agreement 
with the inviscid theory and were approximately 50 percent below the theory for 
inviscid plus viscous interaction in contrast to the present results. Although 
the ratio of wall enthalpy to stagnation enthalpy of reference 26 was lower 
than that in the present investigation, it is felt that the primary cause for 
the differences in the two sets of data stems from Burke’s Knudsen number being 
several times larger than in the present results. 

Only the inviscid plus strong viscous effects (eq. ( 9 ) and eq. (ll)) were 
correlated in terms of ^ x and plotted in the lower half of figures 6, 

7, and 8. The slopes n were determined and were used to evaluate K^. The 
simple addition of equation ( 9 ) and equation (ll) indicates good agreement with 
the correlated experimental pressure results, but the result was considerably 
lower than the experimental heat -transfer data. 

At the present, there are several theoretical approaches to cover the real 
case where both viscous and inviscid effects occur. One approach is to add 
linearly the two effects. (See refs. 21, and 4l to 4-3*) This linear-addition 
approach should be considered nearly valid only if one effect is small compared 
with the other and even then it should not be considered completely valid. 


Sharp Flat-Plate Boundary-Layer Thickness 

A check of hypersonic similarity based on measured boundary-layer thickness 
is difficult because the definition of boundary-layer thickness is arbitrary. 
However, it was felt that some insight into the sharp -leading-edge flat plate 
could be gained from the present Investigation by measuring the apparent 
boundary-layer thickness from schlieren photographs. 
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Figure 9 shows the boundary -layer thickness plotted against the distance 
from the leading edge. The boundary -layer measurements were made from the 
plate surface to the outer edge of the boundary layer. (See sketch (A).) The 
dark line in figure 4(a) denotes the region of large density gradient normal to 
the surface and corresponds approximately with the surface of the thermal 
boundary layer. The boundary layer shown in figure 4(a) is not sharply defined, 
and, therefore, the measurement is only an approximation. It is felt, however, 
that the data are to a large extent representative of the boundary -layer thick- 
ness. The result of the hypersonic boundary- layer theory expressed by equa- 
tion 6 shows fair agreement with the measured thickness. It can be seen that 
for the sharp plate the boundary layer is merged with the shock layer near the 
leading edge (given previously in refs. 22 and 23 ). Boundary -layer thickness 
in the present study is about O .65 inch at a distance of 7 inches from the 
leading edge. 

The concept of including local pressure distribution through the use of 
standard hypersonic approximations can be shown by simply multiplying the 
boundary-layer-thickness equation (eq. ( 6 )) by the square root of the local 
pressure ratio. (See fig. 10.) The equation becomes 



The bottom part of figure 10 shows the boundary-layer thickness in terms 
of the correlating parameter given by the zero-pressure-gradient theory 
(eq. ( 6 )) plotted against Reynolds number based on plate-surface displacement. 
The correlating parameter decreases with decreasing Reynolds number. At the top 
of figure 10 , a modification of the boundary-layer-thickness parameter has been 
made as suggested by reference 1 6 to include the variation of local pressure 
conditions. Equation (19) was employed and also shows good agreement when the 
results were plotted against Reynolds number. Figure 10 indicates that the 
experimental results are essentially independent of the distance x, as was 
predicted by theory. 


Shock Shapes From Schlieren Photographs 

Figure 4 consists of representative photographs of the different leading- 
edge bluntnesses tested at nearly constant Reynolds and Mach numbers. The 
shock shape is essentially determined by the mutual interaction of the shock 
wave and boundary layer for the sharp plate near the leading edge (fig. 4(a)). 
This merged interaction occurs for R^ = 11 . 5 at the conditions tested, as has 
been noted previously by others for similar small R^ values. The schlieren 
photograph (fig. 4(a)) for the sharp plate is believed to be somewhat similar 
to those obtained by Nagamatsu (refs. 22 to 24) for free -molecule flow 
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conditions. Nagamatsu indicates that for large Knudsen numbers, very small Re- 
values, and high flow velocities, the flow at the leading edge is nearly free 
molecular. This free -molecular flow is followed by a rarefied flow, and 
finally attains at the end of the slip distance the continuum flow regime with 
a velocity at the wall equal to zero. These conditions are illustrated in 
sketch (B). 



It was found in the present investigation on the sharp plate that the 
merged layer slope increased from the leading edge to some maximum value near 
the end of the slip region. The shock wave and boundary layer are both merged 
and their formation was delayed in this region. 

The shock-wave coordinates for the various bluntnesses measured from the 
photographs (see sketch (A)) are presented in figure 11. Qualitatively, the 
shock shapes are s imi lar to those observed by others for hypersonic -flow and 
adiabatic -wall conditions. (See refs. 25 and k^>.) The shock wave for the 
sharp plate appears to be attached. Only slight detachment and standoff occurs 
for the blunt leading -edge plates. 

First -approximation blast-wave theory does not agree with the experimental 
shock shapes for any of the leading -edge configurations. (See fig. (12).) The 
second approximation, in general, provides a much closer prediction of the 
experimental shapes and an even closer prediction was obtained from the theory 
of Cheng and Pallone (ref. 4). As expected, the agreement between theoretical 
blast-wave shapes and the experimental shock shapes improves with an increase in 
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leading-edge bluntness. It is believed that the outward displacement of the 
shock wave downstream is due to the large boundary-layer displacement effects. 
Somewhat similar results have been noticed previously by other investigators 
(refs. 25 and 2 6 ). 


CONCLUSIONS 


An investigation of high Mach number low-density flow over a flat plate 
with variation in leading-edge bluntness has been made. The investigation 
included effects of leading -edge bluntness and viscous interaction at zero 
angle of attack. Pressure and heat-transfer data as well as schlieren photo- 
graphs were obtained in the Langley hotshot tunnel. The nominal test condi- 
tions were a Mach number of approximately 20, a stagnation temperature of 

^4-000° R, and a Reynolds number per foot of approximately 0A6 x 10^. Correla- 
tions and comparisons of the experimental results with hypersonic flow theory 
have led to the following conclusions: 

1. Available theoretical methods are adequate for prediction of the pres- 
sure distributions on a blunt flat plate; however, the measured heat-transfer 
distributions are always somewhat higher than the theory. 

2. The experimental pressure distribution on the sharp leading -edge plate 
can be predicted by the strong interaction theory where continuum flow exists. 

3. The hypersonic boundary-layer theory was found to account satisfac- 
torily for the variations in boundary-layer thickness due to the variations in 
pressure distributions on the sharp plate. 

k. The theory of Cheng and Pallone predicted the shock shapes better than 
did the first- or second -order blast -wave approximation theory. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Admini s t r at ion, 
Va., May 19, 1965 . 
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TABLE I 

COORDINATES AT GAGE LOCATIONS 

referenced from leading-edge base line; 
y referenced from center linej 


Gage 

Thermocouple 

Pressure 

orifice 


x, in. 

y , in. 

x, in. 

Y , in 

i 

0.25 

1.25 

0.25 

-0.25 

2 

.25 

.25 

.25 

-1.25 

3 

• 50 

• 75 

.90 

-75 

4 

• 75 

• 25 

.75 

-25 

5 

1.00 

• 75 

1.0 

-75 

6 

1.25 

• 25 

1.25 

-.25 

7 

1.50 

• 75 

1.90 

-75 

8 

1.75 

• 25 

1.75 

-.25 

9 

2.00 

.75 

2.00 

-75 

10 

2.25 

• 25 

2.25 

-25 

ll 

2.50 

.75 

2.50 

-.75 

12 

2.75 

• 25 

2.75 

-25 

13 

5-00 

• 75 

3.00 

-75 

14 

3.25 

• 25 

3.25 

-25 

15 

3.50 

• 75 

3.50 

-.75 

16 

3.50 

1.25 

3.50 

-1.25 

17 

3.50 

1.75 

3.50 

-1.75 

18 

3.50 

2.25 

3.75 

-.25 

19 

3-75 

2.5 

4.25 

-25 

20 

4.25 

2-5 

4.75 

-.25 

21 

4.75 

2.5 




Leading -edge 
c onf igur at i on 


Pressure 

(* 

gage 1A 

-) 

Pressure 

(* 

gage IB 

) 

X, in. 

in. 

x, in. 

y , in. 

0.197 

0.12 

0.0689 

0.50 

.170 

.12 

.068 

• 50 

.130 

.12 

.066 

• 50 

.080 

.12 

.066 

• 50 


Not shown in figure 2. 


ro 

o\ 


TABLE II.- CORRECTION FACTORS AS A FUNCTION OF TIME 

[Subscript vel denotes velocity corrected data; 
subscript cal denotes uncorrected dataTj 


Time, 

msec 

( p arc)vel 

(Parc)cal 

( T arc)vel 

(^arc)cal 

( ^00 )vel 

(Poo )cal 

( p oo)vel 
(Poo)cal 

( T oo)vel 
(■^oo)cal 

( V »)vel 

(^oo)cal 

(^oo)cal 

( ^00) vel 
(^oo)cal 

(^t )cal 

0 

10 

1.45 

1.46 

0.66 

.75 

0.95 

1.39 

0.69 

0.85 

1.03 

1.88 

0.74 

20 

1.30 

.77 

• 97 

I.36 

.72 

.86 

1.02 

1.78 

.76 

30 

1.28 

.78 

.94 

1.33 

.71 

.87 

1.03 

1.78 

.78 

40 

1.27 

• -79 

•95 

1.30 

•75 

.88 

1.03 

1.70 

•79 

50 

1.25 

.81 

• 96 

1.26 

.76 

•90 

1.02 

1.59 

.81 

60 

1.22 

.83 

• 97 

, 1.23 

; .79 

• 91 

1.02 

1.51 

.83 

70 

1 1.19 

1 .86 

■ -98 

i 1.19 

! .82 

.92 

1.01 

1.40 

.86 

80 

1.16 

.88 

.98 

1.16 

.85 

• 95 

1.01 

1.33 

.88 

90 

1.13 

.90 

• 99 

1.13 

.87 

.94 

1.00 

1.28 

• 89 

100 

1.12 

• 93 

1.0 

( 1.10 

: .91 

• 95 

1.00 

1.19 

.92 





(a) Model. 


L- 64- 6875 



(b) Underside of instrumented surface. L- 63- 8643 

Figure 1.- Photograph of model and instrumented surface. 
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Alternate leading edges 
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Figure 2.- Flat-plate model. All dimensions are in inches unless otherwise noted. 
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Distance from leading edge, x, in. 




Distance from leading edge, x, in. 


Figure 3.- Typical variation and reproducibility of measured local pressure and heat transfer on sharp flat plate, a = 0°; ip = 0°; tfe = 0.0003 in. 
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(d) Bfunt- nose flat plate; tj e = 0.080 in 
Figure 4.- Concluded. 





Figure 5.- Pressure and heat-transfer distribution on sharp flat plate. M«, = 20; tj e = 0.0003 in. 








Figure 7.- Pressure and heat-transfer distribution on blunt-nosed flat plate. Mco = 20; t/ e = 0.030 in. 




Boundary-layer thickness 
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